Random backaction in tunneling of single electrons through nanostructures 
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We derive an n-resolved Master equation for quantum transport that includes a dependence on 
the number n of tunneled electrons in system parameters such as tunnel rates and energy levels. We 
apply the formalism to describe dynamical changes due to random backaction effects, for example 
due to local fluctuations of the electrostatic landscape during the transport process. We quantify 
the amount of additional noise on top of electron shot noise due to these fluctuations by giving 
explicit expressions both for sequential and coherent tunneling examples. 
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I. INTRODUCTION 

Single electron counting experimentsi have opened a 
new view on electronic transport through nanostructures. 
Fluctuations of the current can provide insight into inter- 
nal processes that contribute to quantum transport, such 
as quantum coherent oscillations between different parts 
of the nanostructure^, interactions with the environment, 
or correlations leading to non-Markovian dynamics^^— . 

Counting of individual charges during a stationary 
transport process is also fascinating from a fundamen- 
tal point of view, as it touches the question of measure- 
ment of a quantum system in contact with large reservoirs 
(source and drain). The standard way to model a count- 
ing device in transport is the inclusion of counting fields 
X that appear in a natural way as Fourier counterparts 
of n, the number of electrons tunneled during a certain 
period of time t. 

Transport Master equations including x can be derived 
in various ways. An early approach stemming from quan- 
tum optics is the counting of individual stochastic jump 
events^, e.g. in quantum trajectories^. Alternatively, one 
can start from the projection postulate for N, a num- 
ber operator in one of the reservoir, and follow its time- 
evolution* which leads to the Levitov-Lesovik-formulaS 
for non-interacting, or ^-dependent Master equations for 
weakly coupled and interacting systems^. Yet another 
possibility, leading to n-resolved Master equations, is 
the use of the total system plus reservoir wave func- 
tion and perturbation theorj*ii, or the inclusion of ad- 
ditional counter variables in the part of the Hamiltonian 
that describes the coupling between the system and the 
reservoirsi^. 

A somewhat non-trivial aspect appears when the quan- 
tities that determine the transport, such as tunnel rates 
or energy levels, acquire a dependence on n and thus 
become dynamical quantities themselves. In this paper, 
we show that such a dependence can in fact be obtained 
from a model Hamiltonian that includes appropriate mi- 
croscopic interaction terms depending on N ^ the num- 
ber operator for one of the reservoirs that acts as a de- 
tector. By applying the Nakajima-Zwanzig projection 
method^, we extend the usual derivation of the trans- 
port Master equation and obtain an equation of motion 



for p*^"^ , the reduced system density operator conditioned 
on the charge n in the detector. 

We apply our method to model backaction effects that 
have the form of an un-controlled modification of the elec- 
trostatic potential landscape during the transport pro- 
cess. For example, we can account for saturation effects 
due to a local piling-up of charge leading to non-constant 
tunnel rates r„ that change with n. Rather than calcu- 
lating particular model forms of r„, however, we con- 
sider small, but random fluctuations of tunnel barriers 
and confinement potentials in the course of electrons tun- 
neling through the nanostructure and thus assume ran- 
dom, n-dependent parameters in the n-resolved Master 
equation. 

Some generic properties of eigenvalue distributions of 
random transition-rates matrices for rate equations have 
been studied recentlyli. Our focus here is on transport 
quantities such as the stationary current and the Fano 
factor in cases of Master equation parameters that ran- 
domly change with every tunneled particle. We show that 
the additional noise introduced by this randomness can 
lead to significant modifications of current cumulants, in 
particular in the coherent tunneling regime. 

Our formalism is also relevant for the description of 
feedback and control operations in the form of back- 
action from the detector onto the system during the 
transport process, whereby the detector is upgraded 
into an actuator—. In feedback control, the Wiseman- 
Milburn scheme^^ offers a simple way to include in- 
stantaneous feedback in the form of properly re-defined 
jump superoperators in the usual transport Master equa- 
tions. Other ways to include feedback are possible, how- 
ever. Our derivations below provide a microscopic back- 
ground to the recent combination of open loop control, 
i.e. a time-dependence in system and reservoir parame- 
ters, with the n-dependent output of a counting device 
to form a feedback loop^^. 

The outline of this paper is as follows: section II in- 
troduces the n-resolved transport Master equation via 
the Nakajima-Zwanzig method. Section III derives the 
general expression for the ensemble average of the elec- 
tron current and the Fano factor, and in Section IV we 
discuss applications to sequential and coherent tunneling 
through quantum dots. 
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II. MODEL AND DERIVATION OF MASTER 
EQUATION 

Let us consider a situation where electrons are trans- 
feree! from a source ('left') reservoir through a system 
into a drain ('right') reservoir. We use a decomposition 
of the total Hamiltonian into system (S), bath (left and 
right reservoir, B), and systcm-bath interaction (SB), 



H^Hs + Hb + Hsb- 



(1) 



The dynamics of the density operator g of the to- 
tal system (including the reservoirs) is described by the 
Liouville-von Neumann Equation dtQ — —i\H,Q\. We 
trace out the reservoir degrees of freedom, keeping only 
the system and the number n of particle transfered from 
the left into the right reservoir as relevant degrees of free- 
dom. This can be achieved with the method of Nakajima 
and Zwanzig^^ by defining a projector P 

pg = ^p(n)^^in)^^ (2) 
n 

which projects onto the subspace of relevant information. 
Here, 



,(n)=Tr(5'^)(^) = XrB(P„f?F„) 



(3) 



is the reduced system density matrix, conditioned upon 
the number n of electron transfers. Furthermore, the 
trace operation Tr^ ^ (..) over the bath includes the pro- 
jection operator P„ that projects onto the subspace of n 
electrons in the Hilbert space of the right reservoir, and 



(n ) _ PnQBflPn 

0B,a 



(4) 



TTBiPngBfiPn)' 

where g^.o is a fixed equilibrium density matrix of the 
bath. Using the standard Nakajima- Zwanzig technique 
with the projection operator P, we obtain an equation of 
motion for the projected density matrix to second order 
in HsB 



dtPg{t) = PCoPg{t) 



e 



dt' Tr 



(n) 



Pg{t),nsB{t-t') 



(5) 

(n) 



with the Liouvillian Cq corresponding to system and bath 
without interactions. 

To make further progress, we need to specify the 
system- bath interaction Hsb- For the following deriva- 
tion, we assume a simple Hamiltonian for tunneling be- 
tween a quantum dot attached to lead reservoirs. 



nsB - cL^k»m{N)d^ + H.C., 



(6) 



is the reservoir Hamiltonian, and dm annihilates an elec- 
tron in system (dot) level m. In the tunnel matrix ele- 
ment Vkam{N) (defined as a Taylor series in iV), we al- 
low a dependence on the number operator N of the right 
reservoir. This dependence describes the 'backaction' of 
the right detector reservoir onto the tunnel process and 
will render the tunnel rates below to be n-dependent 
quantities. 

In the following, for simplicity we assume a single sys- 
tem level only and thus suppress the index m. The 
multiple-level case can be done in a completely analogous 
fashion. Care has to be taken when evaluating the dou- 
ble commutators in Eq. ([S]). We assume the infinite bias 
limit and thus a bath density matrix gsfi correspond- 
ing to left and right chemical potentials ^ ^ ^r. 
The derivation of the Master equation in the usual Born- 
Markov approximation is now straightforward algebra, 
cf. Appendix A. Setting the single level energy to zero, 
we arrive at 



- (t) , ddi' } - r^"^ d 



(8) 



where the n-dependent tunnel rates are given by the ex- 
plicit expressions 



p(ri) 



27r^|yfe„(n)|2,5(e-£fe„), 



(9) 



and where Ska are the single particle energies in the reser- 
voirs a. We note that Eq. ([8]) generalises the usual Mas- 
ter equation for the single resonant level model which 
is known to be a reliable description of single electron 
transport^- in the infinite bias limit employed here. The 
above derivation can now be easily extended to cases 
where not only the tunnel Hamiltonian Hsb but also the 
system Hamiltonian T-Ls has parameters that depend on 
the particle number operator N of the counter reservoir. 
The n-resolved Master equation then has the generic 
form 



(10) 



where we already used the decomposition of the Liouvil- 
lian into the jump-superoperator that describes 
tunneling of electrons out of the system, and the non- 
jump part >Cg"' . Note that the appearance of indices 
n and n ~ 1 reflects the fact that we are within the 
Markovian approximation (there are no indices n — 2 etc.) 
and only consider unidirectional transport (there are no 
'back' jump terms with index n + I). 



where cJ,q, is the creation operator for an electron in sin- 
gle particle state k with energy Eka in lead a =L/R 
(left /right). 



(7) 



k,a=L,R 



III. RANDOM ENSEMBLES 



The n-resolved Master equation Eq. (jlOp can in prin- 
ciple be solved easily on a computer once the matrix ele- 
ments of the super-operators £q"'' and are known. 
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The easiest situation is, of course, the case where the 
super-operators do not depend on n at all, and a simple 
Fourier transformation of Eq. ([T0|) according to 



where we defined the propagator 



Px 



(11) 



simply leads to a single, n-independent equation for 



(12) 



In the following, we assume that with each single elec- 
tron tunneling event, the complete configuration of the 
total system changes in a stochastic manner. For exam- 
ple, charges locally piling up near one the tunnel junc- 
tions will modify the electrostatic potential of the tunnel 
barriers and also affect the confinement potential of the 
nanostructure. In general, these effects will lead to a 
dynamical modification of tunnel rates or energy levels, 
depending on n, the number of charges tunneled through 
the structure. Usually, this n-dependence is small and 
assumed to be negligeable, and the electrostatics of the 
total system is described by means of some average po- 
tential landscape thereby giving rise to constant rates 
and energies. 

A more realistic approach is instead to allow for an 
n-dependent variation of the parameters in the Master 
equation. As it is practically impossible to microscopi- 
cally model the dependence on n, the assumption of a 
random variation of these parameters should be a good 
starting point. It is therefore physically reasonable that 
the n-dependent super-operators in Eq. ([TU)) are indepen- 
dent, statistically distributed and uncorrelated to each 
other. The total Hamiltonian H underlying Eq. ((TU]) 
thereby becomes random, as the Master equation and 
the density matrix [t) itself. Expectation values cal- 
culated with the help of p'-"-' {t) then involve two averages, 
i.e. the usual trace average (for a given realization within 
the random ensemble), and the 'disorder' (ensemble) av- 
erage that we denote by (. . .) in the following. 

The probability p{n, t) of n electrons having tunneled 
into the right reservoir until time t (counting is started 
at t = 0) is given by 



p(n,t)=Tr(p(")(t)). 



(13) 



As this quantity is in general quite difficult to calculate, 
we will be satisfied with a calculation of mean values and 
variances in most of what follows. 



A. Ensemble Average 

To calculate the ensemble average of p^"^(t), we first 
Laplace transform and recursively solve Eq. (jlOp . which 
leads to 

p„(z)=7'(")(z)j("'i)...J«p(i)(z)j(")7^(")(z)po, 

(14) 



(15) 



and po is an initial condition. In the long-time limit 
t -> cxo for the stationary state discussed below, the 
choice of po becomes irrelevant in our Markovian theory. 

We now assume the distribution function of the n- 
dependent super-operators in Eq. (jlOp to factorize into 
independent and identical distributions for all n. This 
is a reasonable assumption, as long as strong inhomo- 
geneities are excluded. A counterexample where this as- 
sumption fails would be saturation effects, e.g. a gradual 
and steady increase (or decrease) of tunnel rates with 
increasing number n of tunneled electrons. 

Using the factorization assumption, the expectation 
value for the system density matrix in Laplace space 
reads 

(p„(z)) = W(z)"P(z)po, (16) 
where we defined the ensemble-averaged super-operators 
W(z) = (j("¥(")(z)), P(z) = (7^(")(z)). (17) 

Here, the product JV in the definition of Wlz) is anal- 
ogous to the product of superoperators used previously^ 
in a discussion of waiting time distributions for Marko- 
vian quantum transport. 

Using the Fourier transformation Eq. ljlip with respect 
to the counting field Xi 'we obtain the ensemble average 



(Px(z))= l-e^^W{z) P{z)po 



(18) 



From this representation one can now directly obtain the 
ensemble averaged moments of the Full Counting Statis- 
tics in Laplace space, defined as 

{mk{z))^i-^TTd'^{Mz))\^^^, (19) 
where we set the electron charge to unity. 



B. Current and Fano Factor 

The long-time limit for the ensemble averaged electron 
current is defined as 



(loo) = ( lim dtmiit)). 

t^OQ 



(20) 



This can be obtained directly in Laplace space via 

(/oo) = limz2(mi(z)), (21) 

where here and in the following the hat denotes that 
the quantity was Laplace transformed, i.e. f{z) — 
dte^^*/(t), and the correspondence between the limit 
lim f{t) in time domain and lim zf{z) in Laplace space 

t— >oo z— >0 

was used. 
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Some more effort is needed to obtain the long-time 
value of the Fano factor 



(Foo) = lim 



{m2{t)) - (mi(t))- 

(TOl(i)) 



(22) 



Here, a direct calculation in Laplace space can not be 
efficiently executed. Fortunately, the Laplace transfor- 
mation has the convenient property that terms of the 
order z"'' in Laplace space correspond to terms in 
the time domain. We thus derive the Fano factor by first 
performing an expansion of m2(z) in inverse powers of 
z up to order = 1, after which we apply the inverse 
Laplace transformation in order to extract the long-time 
limit. 

This procedure can also be applied in an efficient way 
for the higher cumulants^-. In the following, however, 
we restrict ourselves to the first two cumulants, i.e. the 
average current and the Fano factor, which we evaluate 
and discuss for some representative and simple quantum 
systems. 



IV. EXAMPLES 

We now evaluate our results in detail for single electron 
tunneling in representative examples of sequential and 
coherent tunneling. 



We start with the calculation of the first moment, 
which evaluates to (toi(z)) = — \, and the long- 



time limit of the current is therefore 



(/oo) = 



^(0) 



(26) 



This result is interesting in that it yields the current as 
the inverse of the ensemble averaged time r and not, as 
one might first have expected, the ensemble average (F). 
In the 'clean' case, i.e. without random variations and a 
delta function distribution of the rates, /(F') = (5(r' — F), 
the current is indeed given by (/oo) = Tj and t = F^^ is 
the mean value of the exponential distribution of waiting 
timcs^^ for the Poissonian process describing the tunnel- 
ing. As we are dealing with single electron tunneling, the 
waiting timei^i^i"— , rather the tunnel rate, is in fact the 
more fundamental quantity to describe the average cur- 
rent: if it takes on average a time r to send a single object 
from place A to place B, on average a stationary current 
of 1/r objects flows between A and B. This definition 
also holds if there are random variations of the time t in 
which case t has to be replaced with its ensemble average 

Next, we calculate the Fano factor, for which we first 
need the ensemble averaged second moment 



TO2(z) 



z^VizY z^V{z) 



(27) 



A. Tunnel Contact 

The simplest case is a tunnel contact in which the 
left (source) reservoir is connected to the right (counter) 
reservoir without anything in between, i.e Hs = 0, which 



corresponds to super-operators £g 



(n) 



defined by scalar rates F^"' in our Master equa- 
tion Eq. (flUl) . The ensemble averages in Eq. (IT71) then 
become 



Viz) 



1 



z + T 



VV{z) = 1 - zV{z), 



(23) 



We insert the Taylor-expansion of Viz) around z = 0, 
Viz) « (r) - z(t2) + z'^{t^) + 0(z3), into 1^ and the 
expression for mi(z). After Laplace back-transforming, 
the moments are then combined to give the Fano factor 
in the long-time limit. 



(r)2 



-1 = 1 



. Var(r) 



(28) 



where Var(r) = (t^) — (r)'^ > is the variance of t 
which thus renders the ensemble averaged Fano factor as 
greater than unity. 



where F denotes the random tunnel rate of the junction. 
Note that 'P(z) is the generating function for moments 
of the inverse tunnel rate t = F~^, 



^.^fc+ij, ^ (^_^^k9z'P{z) 



(fc) 



(24) 



Using Eq. (fTS)) and the initial condition po = 1 (fixed by 
the normalization) , we obtain 



V{z) 



1 + zViz) - 1 



(zV{z 



(25) 



B. Single Quantum Dot and Ring 

The tunnel junction example above can be easily gen- 
eralized by formally combining K > \ junctions in series 
but still only allowing for one electron tunneling through 
the system at a time. The case K = 2, for example, 
corresponds to a single- level quantum dot, a system that 
has been well studied experimentally in the pasljii^. In 
general, a sequence of K tunneling contacts can be re- 
garded as a ringiS when the transitions between the first 

(n) 

K states, 1^2---^K^\aX rates F^ ' are in- 
terpreted as no-jump superoperators and the transition 
-ftT as the jump process. The n-dependent Liouvillian 
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super-operators are 



-^0 ^ 



/ _p(«) 








p(n) 
^ 1 


^ 2 








p(n) 
^ 2 


p(n) 
^ 3 



V 



.(") 



.(") 

K 



J 



(n-1) 



(29) 



with {{K\ = (0, and = (1,0,...)^- 

Correspondingly, the ensemble averaged super- 
operators W{z) and V{z) now become K x if-matrices, 



(Wni Wn2 Wm 




W(z) 



Wyik \ 








(30) 








0/ 




Pi 











WIP2 


P2 








WIW2P3 


W2P3 


^3 





V 



Wk 



Wji2Pk 

Wk 



V{z) 



, . . Wk^iPk Pk ) 

with scalar parameters Wi — = (l^T") ^^'^ 

Wm = Ylf^i ■ The disorder averaged moment gener- 
ating function now has to be defined as usual via the trace 
operation, which for an initially empty system yields 



K i 

i=i k=i 

1 - e'xWni 



(31) 



After back-transformation into the time domain and tak- 
ing the long-time limit, current and Fano factor read 



(^00) 



1 



K 

i=l 



i=l 



2Var(T,)) 



K 

i=l 



(32) 



Again, the result for the stationary current has its sim- 
ple origin in an average sequential passage time through 
the system, for example {loo) — {tl)+{tr} ^^^'^ single 
level dot {K = 2). The ensemble averaged Fano factor in 
this case reads 



{tl)^ + {tb)^ , 2Var(rL) + 2Var(rfl) 



{{rn) + {tlW 



(33) 



This result has a very natural interpretation: the first 
term in (i^oo) is identical with the 'clean' (non-random) 
Fano factor for a single level dot with the inverse rates 



and F^^ replaced by the ensemble averaged single- 
junction waiting times. The second term describes the 
increase of the current fluctuations that are due to the 
fiuctuations of the waiting times t^j and in the random 
ensemble. 

We also note that the Fano factor in the sequential 
tunneling case obeys the inequality 



{Foo) > 



i=l 
K 



(34) 



where the last inequality simply follows from the Cauchy- 
Schwarz inequality {er)'^ < ||e|P||'r|p = if ||t|P, where r 
is the vector with the K components (t)^ and e the vector 
with the K components 1. 



C. Double Quantum Dot 

In contrast to the sequential tunneling examples of 
the previous section, a coherently coupled system of 
two quantum dots (double quantum dot, DQD) can- 
not be described as a simple serial combination of sin- 
gle dotaiii^"— . The DQD Hamiltonian describes a sin- 
gle additional, spin-polarized electron on two levels \L) 
(left dot) and \R) (right dot) that have energy difference 
e and are coherently coupled by a matrix element Tc- 
Using pseudo-spin Pauh matrices (Tz = \L){L\ — \R){R\, 
ax = \L){R\ + \R){L\, the Hamiltonian reads 



n = Hs+H,es+HT 



(35) 



= (cL^fco(iV)|0)(a|-h i?.c.) 



(36) 



k,a=L,R 



with the 'empty' state |0) and the standard tunnel Hamil- 
tonian T-Lt for coupling to the left and right reservoirs 
Hres- Here, we assume a dependence of the system part 
Tc on the counting number operator N in addition to 
this dependence in Vka{N). 

The n-dependent superoperators of the DQD in the 
basis p — (po, PL, PR, ^PRL, '^Prl) then have the form 



-^0 ~ 





^) 











-n(" 
^ L 



















^ R 





-2T(^" 











2 


— e 


V 


rp(n) 


rp{n) 


e 


^ R 

2 


^ R. 


R)){{R\ 









2 / 

(37) 

with the non-jump superoperator >Cg"^ and the jump su- 
peroperator with ((^| = (0,0,1,0,0) and \R)} = 
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(1,0,0,0,0)'^. Furthermore, T^"^ is the eigenvalue of 
Tc{N) corresponding to n electrons in the right lead, 
and r^"'' , r^'"* are n-dependent tunnel rates defined as 
before in the sequential tunneling case, Eq.®. 

The ensemble averaged stationary current is obtained 
from the moment generating function (cf. Appendix B, 
Eq. (IBip ) and evaluates to 



1 



{rL) + {rn) (e2(r|) + 2) + i(r«)(r2) ^ 
where we defined 



(38) 



(ri;) = (^).(39) 



Some interesting observations can be made from the ex- 
plicit expressions for (loo)- First, Eq. (|38l) reduces to 
the known resul1i^iS£iS^ in the clean, non-random limit 
for the parameters Tl, ^r, and Tc- 

Second, we recognize that in the expression for (loo), 
the random coupling to the right reservoir enters in the 
form of two independent averages, i.e. the mean values 
(Tfi) and (tr) of the right tunnel rate and its inverse, re- 
spectively. In other words, when fixing the averages Eq. 
([M]) there is in general no co- incidence between (loo) and 
the corresponding current loo without random fiucuta- 
tions in the parameters. For a random ensemble of right 
tunnel rates Tn, the ensemble averaged current (loo) is 
always smaller than the corresponding clean result /qo , 



1 



< 1, 



{rn){Tj,)-l. (40) 



Here, the parameter a depends on the probability distri- 
bution, but it is always non-negative owing to Jensen's 
inequality g{{x)) < {g{x)) for convex functions g{x) in 
the special case g{x) — i. 

For example, if we consider a uniform probability den- 
sity distribution 



i -§ + {rR)<x<§ + {Tii) 
else 



(41) 



with average (tr) and width a for the inverse right tun- 

nel rate tr = I/Tr, we find {Tr) = i log (l - ^^^) ■ 

For a — > 2{Tfj), i.e. when arbitrary small right waiting 
times tr — >■ and thus infinitely large tunnel rates F^; 
become possible, the current (loo), Eq. (|38|) . becomes 
more and more suppressed, cf. Fig. [1] This suppres- 
sion of (/oo) is a manifestation of the Zeno effectj^l, i.e. 
strong detection via the right counter reservoir though at 
fiuctuating detection strengths. It also occurs for other 
forms of probability density distributions /^ (^^.^^{x) that 
give sufficient weight to small waiting times tr. 

Similar to the ensemble averaged current, we also ob- 
tain the Fano factor from the moment generating func- 
tion, cf. Eq. (jB4[) . where second moments with both (r|.) 
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FIG. 1: Stationary currents (/oo) (lower part, Eq. ((38])), and 
Fano factors {Foo) (upper part, Eq. (|B4|) ) as a function of 
internal dot bias e for non-random case {a = 0) and random 
distributions with width parameter a > (inset, Eq. (|4ip ') 
of inverse tunnel rates tr ~ at right barrier of double 
quantum dot. Average (tr) — 5, other parameters fixed: 
internal bias e, tunnel coupling Tc ~ ^ (in units of left tunnel 
rate Tl)- 



and (F|j) enter (in analogy to the parameter a, it is pos- 
sible to introduce positive parameters which describe the 
relation between moments (t*^) and (F*^) for k > 1—). In 
contrast to the sequential tunneling result Eq. p3p . one 
can no longer separate the fiuctuations of the current due 
to the stochastic tunneling process and due to the fluc- 
tuations in the tunnel rates, which is due to the coherent 
coupling between the two dots. Such a separation is only 
possible at very large internal bias e, where 



lim (Foo) 



+ 2 



, o var(rfl) 

+ {tr}' 

var(T.^) var(rfl)var(T|,) 



{rR)H4y 



(42) 



In this limit, the transport is essentially determined by 
the single Poissonian process at the right barrier and 
we recover the Fano factor from the sequential tunnel- 
ing case, Eq. (1^ . plus the additional contribution due 
to the fiuctuations of as described by var(T|,). 

For smaller e, the Fano factor as a function of e has 
no longer a line shape that is simply shifted as compared 
with the clean (non-random) case, cf. Fig. [T] The box 
distribution example, Eq. (|411) . also shows that (Foo) 
strongly grows at e = when small tr become likely and 
fluctuations between very long and shorter waiting times 
become stronger. 
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V. DISCUSSION AND CONCLUSION 

For sequential tunneling, the waiting times and their 
random distributions determine the ensemble averaged 
transport quantities, which is shown already in the sim- 
plest example of a single tunnel contact where the current 
is given by the inverse ensemble average {t)~^ and not 
the averaged rate F = r~^, cf. Eg. (|26| . The random 
fluctuations of the across the barriers simply give an 
additional contribution on top of the usual shot noise 
calculated with ensemble averaged waiting times, cf. Eq. 
([32|) . We expect this correction to be small for all current 
cumulants, as long as the random fluctuations are small, 
and experimental results^^ indeed show that a modeling 
with fixed Master equation parameters is very successful. 

For coherent tunneling, however, the situation seems to 
be more involved. There can be additional n-dependent 
fluctuations of internal system parameters (energies, tun- 
nel couplings), and we find that fluctuations of the wait- 
ing times and the inverse waiting times (tunnel rates) de- 
termine the transport, cf. Eq. (I55|) for the DQD current. 
We numerically confirmed that the modifications of DQD 
current and Fano factor are small if the waiting times tji 
at the right barrier fluctuate only weakly. In contrast to 
the sequential tunneling case, however, the noise due to 
disorder in the parameters is not simply additive. In par- 
ticular for higher cumulants, this might affect a straight- 
forward interpretation of experimental data based on a 
modeling with Master equations that neglect the kind of 



J 



random backaction discussed here. In contrast, our for- 
malism allows one to take into account these effects, how- 
ever at the expense of introducing at least one additional 
parameter characterizing a distribution function such as 
the width a in Eq. (|4T|) . Quantitative statements (such 
as Eq. (|38|) or Eq. ()B4|) ) then require the Liouvillian of 
the particular nanostructure under consideration. 

A further application of our method would be the cal- 
culation of noise spectra. An interesting line could be the 
modeling of I// charge noise contributings to transport 
with dynamically changing tunnel rates F„. 

Finally, another open (though quite challenging) as- 
pect is to go beyond the Markovian scheme used in this 
paper. One could then test in how far random parame- 
ter fluctuations effectively would wash out the quantum 
memory effects in non-Mar kovian noise features^. 

This work was supported by DFG grant BR 1528/7-1, 
1528/8-1, SEE 910, GRK 1558, the Heraeus foundation, 
and the DAAD. Discussions with W. Belzig are acknowl- 
edged. 



Appendix A: n-resolved Projection Method for 
Master Equation 

The derivation of the Master equation starts from the 
double commutator expression Eq. ([5]) and the form Eq. 
dH]) for the system-bath interaction Hsb- Specifying to a 
single dot level for simplicity, we obtain the eight terms 



= Tr„(") 



©p' 



ka 



(Al) 



\kak' a' n' / \kak' a' n' 



,J2cl^VUN)d + H.c j 

ka. J / 

'^dS+TrB^''^( E ^^s>/.^W4acL'^4'o'(^))p'"'^rf^d 

tk'a'n' / 

E VL{N)d,Js2cl^>V,,^,{N)\ d^p(-">d 

tk' a' 7i' / 



- Trs^-M E 4a^^-(^)^'B>;'a'(^Kw)rfV"'^rf^-TrB(") 

\kak' a' n' / \kak' a' 7i' 

/ \kak' a'n' / 

E ^/r'o'WCfc'a'4a^^aW^?S)rf^'^V"'^+TrB(")[ E 4'a'H'oKiV)l^;jiV)5,„4";A dJt^(n') 
\kak' a' n' / \kak' a'n' / 



- Trs^") E Vk'^'iNK'a'B^B 

\kak' a' n' 



Trs 



where the tilde in the operators abbreviates the time- 
dependence, e.g. Cka = CfcaC"*^''"*^*"* where tka is a 
single-particle energy in lead a, and p'" ^ depends on t' . 
We assume the infinite bias limit and thus a bath den- 
sity matrix qb.{) corresponding to left and right chemical 



potentials pL ^ Q ^ fJ-R- The bath correlation function 
for the jump term for the right (counter) reservoir, e.g.. 
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is then determined by 



(A2) 



uncorrelated distributions for F^"'', F^'' and Tq"'' that 
are independent for ah n. The calculation of the mo- 
ment generation function is analogous to the calculation 
for the sequential case, with the result 



where we kept in mind that P„ projects onto the sub- 
space of n electrons in the right (and not the left) reser- 
voir a = i?, and N in Vk'a'iN) is the number opera- 
tor of the right reservoir. The expression Eq. (jA2[) is 
still approximate in the sense that it neglects the differ- 
ence between the usual grand-canonical Fermi functions 
fkR = Tre QB,oc\piCkR and the canonical projections of 
the occupations at fixed particle number n. This differ- 
ence, however, becomes irrelevant for fiR — > — oo, or more 
generally, if for arbitrary chemical potential ^lr the pro- 
jection P„ is on a subspace with n + Nq particles, where 
A'o is a macroscopic particle number determined by ^r, 
in which case the a = R term in Eq. (|A2[) would have to 
be replaced by the Pauli block factor 1 — fkR- 

A non-jump-term associated with the left (source) 
reservoir (again for /i^ 3> 3> (J-r) has the form 



TrB (Pn Q^B.l cL' 14'a' {N)V^c {N)Ck^ Pn ) 
= 5nn'5kk'\Vka{n)\'^5a,L, 



(A3) 



M{ix) 



Tr I {[z - ^0 



{A{z))-Wl{z){B{z))-\Pl{z) 
exWL{z){TRB{z)) 



Po 



(Bl) 



where again an initially empty dot was assumed and 



1 



-), Wl = 1~zPl{z) (B2) 



z + Tl , 

refer to the left tunnel barrier and are defined in analogy 
to Eq. (1^ . whereas 



where we used the fact that N commutes with c^^l . Cru- 
cially, the tunneling matrix elements now depend on the 
eigenvalues n of the number operator TV. In the infinite 
bias limit, the Markov approximation becomes exact in 
the integration over time t' , and the time dependence in 
the Fermion operators leads to a delta function that to- 
gether with the sum over k yields the n-dcpendcnt tunnel 
Eq. (HI). 



rates F^" 



A{z) = 



{Tr + 2z) {{Tr + z)iTR + 2z) + ATI) + 4e'(Ffl + z) 



B{z) ^Tt 



{Tr + 2z)2 {z(Tr + z) + 4T2) + 4ze2(Ffl + z) 

Tr + 1z 

[Tr + 2z)2 {z{Tr + z) + 4r2) + Aze'^iTR -f z) 

(B3) 



Appendix B: Moment Generation Function and 
Fano Factor for Double Dot 

Using the techniques outlined in section III, we per- 
form a random ensemble average where we again assume 



are combinations containing the random variables T r 
and Tc that have to be averaged over in the definition 
Eq. (|B1[) . Using these result, we obtain the stationary, 
ensemble averaged Fano factor for the double quantum 
dot 



{loo? 



{tlY + 2{Ti) + (r^) ( -Ae^TRy + IQe^r'^) - {Tr){tr) + - 



(B4) 



{r^T? [^'{rRy + \^'{^R){rR) + YQ{TR?)+{4) [^^\rl) + \{TI) + e'^ - ^tr^ 



rl) 
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